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Gauge invariant density-density correlators yield detailed information on hadron structure. 
Hadron deformation and form factors can be extracted for momentum transfers up to about 
6 GeV 2 . We use stochastic techniques and dilution to compute the all to all propagator required 
for the exact evaluation of density-density correlators. We present first results for the pion form 
factor. 
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1. Introduction 

A number of experiments have been performed recently to determine a possible deformation 
in the nucleon [jl|]. Lattice QCD provides a model independent method for studying hadron defor- 
mation via the evaluation of density-density correlators. We show that it is now feasible to evaluate 
density-density correlators for mesons and baryons without any approximation. Such an evaluation 
requires a computation of all to all propagators, which is carried out using stochastic techniques 
and dilution. We show results for hadron charge distributions and demonstrate how hadron form 
factors can be extracted from density-density correlators. We present first results for the pion form 
factor. 



2. Density-density correlators 

For general time insertions the density-density correlator is given by: 

C(y,h,ti) = J d 3 xdh(h(z,t)\f {x + y,t 2 )M(xA)\h(0,0)^ (2.1) 

where j is the time ordered density operator i.e. yAx) =: qf(x)joqf(x) : for a quark of flavor /, 
while \h) denotes any hadronic state. The integration over the x coordinate sets the sink momentum 
equal to that of the source. The integration over z sets the momentum of the source and sink to 
zero. These two sums require summation over both spatial coordinates of the quark propagator 
connecting the density insertions with the sink and thus one needs to evaluate all to all propagators. 

In this work we use stochastic and dilution techniques to compute all to all propagators. It has 
been demonstrated [§] that an estimate for the all to all propagator can be computed by inverting 
the Dirac operator for an ensemble of noise vectors as sources. If these noise sources are created 
in such a way that they obey: 

(r}l{x)r} b v \y))=8(x-y)8 c , h 8^ v and {^{x)) r = (2.2) 

where r denotes the number of noise vectors in the ensemble, then one can invert for each of these 
noise vectors and obtain an estimate for the all to all propagator: 

(Y h v (xK l (y)) r ^(M-\x,y)) b v ^ (2.3) 

where y/^(x) = (M _1 (x^^v^rj^y) is the individual solution for each of the r noise vectors. We 
additionally employ a method known as dilution, which gives a better estimate for the all to all 
propagator The method essentially proposes a way for creating the noise source ensemble in 



Eq. ([2.2D. In this work we use even - odd, color - spin diluted noise vectors i.e. each of the r noise 
vectors of the ensemble has random entries only on either odd or even spatial sites and on one color 
- spin component with all other entries set to zero. Since the time slices of the source, sink and 
density insertions are fixed the noise vectors are in effect already diluted in time. 



3. Wave functions 



We first consider the case where the two density insertions of Eq. ( |Q| ) are taken at equal times. 
This equal time correlator, shown schematically for mesons in Fig. [TJ, reduces to the wave function 
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squared in the non - relativistic limit. For baryons, the wave function is generally a function of two 
relative coordinates and thus requires the calculation of a three-density correlator depicted in Fig. ||. 
In this work we compute the one particle correlator with two density insertions on two different 
quark lines of the baryon. This is equivalent with integrating the three-density correlator over one 
relative coordinate as demonstrated in Ref. |Q] when the diagram in Fig. |2|a was considered. 




p d (X+y,n) (a) (b) 

Figure 1: The density-density Figure 2: The three-density correlator for a baryon. The density 
correlator for a meson. insertions are denoted by the crosses. 

In this study we use 200 unquenched gauge configurations produced by the SESAM collabo- 
ration [||] on a lattice of size 16 3 x 32, at j8 = 5.6 and K = 0.157. Using the mass of the nucleon 
in the chiral limit we obtain aT y = 2.56 GeV yielding m n ~ 950 GeV. We use Wuppertal smearing 
to construct the interpolating fields for the source and sink. HYP smearing is applied to the gauge 
links that enter the Wuppertal smearing function. This decreases considerably the time interval 
needed for the suppression of excited states and allows to place the density insertions as early as 
three time slices from the source [|6|]. Although placing the density insertions close to the source 
or sink is not important for the equal time correlators it is crucial for the extraction of form factors 
discussed in the next section (see discussion connected with Fig. ||). Similarly to satisfy the re- 
quirements for the extraction of form factors, the minimal separation between source and sink that 
is required is t = 14. We take the same source-sink time separation also in the evaluation of equal 
time correlators so that the same set of noise vectors can be used. The evaluation of the density- 
density correlator requires a forward propagator from the source and two all to all propagators, 
namely one from the density insertion at t\ = 3 and one from the sink. We note that the time slice 
of the second density insertion can be varied without additional cost a crucial observation for the 
evaluation of form factors from density correlators. For each stochastic inversion we use 6 sets of 
even - odd, color - spin diluted noise vectors. For comparison a separate computation was carried 
out where no summation was performed on the sink and thus the momentum of the hadronic state 
is not explicitly set to zero. 

3.1 Results 

In Fig. [| we show the density-density correlators for the pion and the p meson. The distribu- 
tions are compared to the case where no explicit zero momentum projection is carried out. If the 
non-relativistic limit is a good enough approximation so that factorization of the center of mass 
momentum can be assumed then there should be no difference between the two evaluations. For 
the pion there is very little difference showing that either factorization is a good approximation or 
exponential suppression of the time evolution sufficiently damps out higher momentum states. On 
the other hand, the p distribution is broader when explicit projection to the state of zero momen- 
tum is carried out. This shows that even for this large quark mass, factorization of the center of 
mass momentum does not strictly hold. These distributions are fitted using an exponential Ansatz 
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Figure 3: The correlators for the pion (bottom) 
and the p, in the zero-spin projection, (top). The 
crosses are the results of the exact evaluation and 
the circles are the results when no explicit zero- 
momentum projection is carried out. 



Figure 4: The correlators for the nucleon (top) and 
the A in the ±3/2 spin projection (bottom). The 
notation is the same as that of Fig. ||. 



allowing the extraction of the root mean square (r.m.s.) radii. For the pion we find an r.m.s. radius 
r,-.m.s = 0.197(2) fm and for the p meson r rm v = 0.576(6) fm. The latter is approximately 25% 
larger as compared to what is obtained without explicit zero momentum projection. 

In Fig. [| we make the same comparison for the nucleon and the A correlators. Both states are 
broader when all to all propagators are used to project to zero momentum. An Ansatz that describes 
these distributions accurately is of the form e^^l^" . Fits to the distributions yield d ~ 1.2 for the 
nucleon and d ~ 1.5 for the A. The r.m.s. radius for the nucleon is only slightly bigger when the 
zero momentum projection is implemented whereas for the A we find a 30% increase. 
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Figure 5: The p (left), nucleon (center) and A (right) correlators projected along the axes. Top: with zero 
momentum projection. Bottom: with no zero momentum projection. 

To probe hadron deformation we project our results for the density-density correlators along 
the spin axis (taken as the z axis) and perpendicular to the spin axis. In Fig. || we show the projec- 
tions for the p meson in the zero spin projection, the nucleon and the A in the ±3/2 spin projection. 
We again compare the exact results with the case where the summation on the sink is omitted. The 
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p correlator shows an elongation along the spin axis. As expected this elongation is much clearer 
when the state is projected to zero momentum. For the nucleon on the other hand, no deformation 
is detected within this method. The A baryon is too noisy to extract any useful information con- 
cerning a possible deformation. More statistics are clearly needed if one is to draw any definite 
conclusion for this hadronic state. 



4. Form factors 

The density-density correlator can be used to extract hadronic form factors [^]. Writing 



Eq. ( |2.1| ) on the lattice and inserting three complete sets of hadronic states, one derives the ex- 



pression: 

|2 



V |(n^(Q)\| 2 IM^ Mjj^MMl^-MM (4.1) 



p.n.s 



where s represents internal degrees of freedom of the hadronic state and J s represents the smeared 
interpolating field used. In Eq. ( |4.1[ ) we assume isospin symmetry (jo = j Q ' = j^). We also assume 
that the time separation between the density insertions from the source and sink are large enough 
so that excited state contributions are suppressed. Now if, in addition, t% — 1\ is large enough and 
we take the Fourier transform over the y coordinate we arrive at the following expression: 

G h ff h ( S ;q; tl ,t 2 ,t) = V I (/, \h s (0) ) I 2 I (°) 1 JO \ h f (fl) f g -s»(g)(fe-rO g -*«Hfc-*)) , (4.2) 

from where the hadron form factor can be extracted in the isospin limit. Note that in practice one 
can calculate the above quantity for any hadronic state with the same set of propagators. Addi- 
tionally one has the form factor for all momentum transfers since these follow from the Fourier 
transform of the density-density correlator. 

To demonstrate that this method works in practice we consider the simple case of the pion. 
The pion as a pseudoscalar meson has only one form factor, which is given by: 

{nip'^Jp = {p„+p'^F % (q 2 ), (4.3) 



where = p' — p^. Combining Eq. (43) and Eq. ( fk2| ) the form factor F K (q ■ ) can be calculated 



by dividing the density-density correlator with an appropriate combination of two-point functions: 

* W ' r-»>l,(£,( e )+M,) 2 G" ?,I 2 -Il )G" {0,1) 

(2-fl»l 

The ratio of two-point functions (G nK ) used in the above equation is the simplest that cancels the 
unknown overlaps of the interpolating fields with the pion state and the exponential time depen- 
dence. However it is not the optimal as far as signal to noise is concerned. Ideally one wants to use 
two-point functions for the shortest possible time separation since the further the source is from the 
sink the noisier the propagator becomes. Thus one must find the optimal values for the insertion 
times and source-sink separation. The dilemma is that although t 2 — 1\ should be large enough to 
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Figure 6: The density-density correlator with the 
insertion time slices kept fixed at ti — t\ = 4. For 
this test we only used point to all propagators. 



Figure 7: The pion form factor for various values 
of the time separation between the density inser- 
tions. We also compare with the prediction using 
vector meson dominance (solid line). 



damp any intermediate excited states, the sink-source time separation t should be as small as pos- 
sible so that the statistical noise on the two-point functions is minimal. Using smearing techniques 
we can ensure that ground state dominance holds for separations of the density insertions from the 
source and sink as small as three time slices. This is demonstrated in Fig. ^ where we plot the 
density-density correlator for the pion with the sink separation fixed at t = 14 and keeping the time 
separation between the insertions at a constant value of t 2 — 1\ =4. We vary the time slice of the 
density insertions by one starting from t\ = 3 and ending at t\ = 7. It can be seen that the results are 
identical for ?2,?i > 3 and therefore we can set the first density insertion at t\ = 3 and the second at 
t2 = 1 1 three time slices from the sink. 

Since t\ = 3 and t — ?2 = 3 are the shortest separations we arrange so that the appropriate ratio 
involves two-point functions in terms of these time separations or time intervals of similar length 
e.g. ^2^- which takes the maximal value of 4. The optimal ratio is then given by: 

F 2 ( 2 )= Um 4£ ff (g)M g G*j j *(q-,t u t 2 ,t)[G™(qA)] 4 G™(0A) 

where G™ and G™ denote smeared - smeared and smeared - local pion two-point functions re- 
spectively. To check for convergence we plot in Fig. ^]f^(^ 2 ), given in Eq. (4.5), as a function of 



?2 — h- As can be seen the pion form factor has converged for t2 — t[ = 8, which is the maximum 
even separation that can be achieved on this lattice. On the same plot we compare the form factor 
to the results obtained assuming vector meson dominance, i.e. taking F Jl (Q 2 ) = X+ Qij m i where 

<2 2 = — q 2 is the Euclidean momentum transfer squared. For m p we take the p meson mass com- 
puted on the lattice which is m p = 1 .270(8) GeV. In Fig. || we compare the pion form factor to other 
recent lattice results extracted using the standard approach of evaluating three-point functions. 

Within this method we obtain reliable results up to momentum transfers of Q 2 ~ 6 Gev 2 . The 
momentum transfers we can extract are limited by the fact that for high momentum transfers the 



density-density correlator given in Eq. (4.2) becomes too noisy, becoming negative and the square 
root cannot be taken. 
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Figure 8: The pion form factor versus Q 2 with similar m K /m p ratios. The triangles and crosses show 
quenched results form Ref. [||], the asterisks are results using twisted mass fermions Q], the squares are 
the results of this work and the circles denote the experimental results. The solid curve is the results 
obtained using vector meson dominance with the physical p meson mass. 



5. Conclusions 

Stochastic techniques combined with dilution are employed in the evaluation of all to all prop- 
agators needed for the exact computation of density-density correlators. The equal time four-point 
correlators reduce to the wave function squared in the non-relativistic limit and thus provide de- 
tailed information on the quark distributions inside the hadrons. In particular the deformation of 
the hadron can be studied. We clearly detect a deformation in the case of the p meson. However, no 
definite conclusion regarding deformation can be reached for the A, at least within these statistics. 
Another application of density-density correlators is the evaluation of hadron form factors. We ap- 
ply this method to compute the pion form factor up to Q 2 ~ 6 GeV 2 . The results we obtain within 
this method are comparable to the results extracted from the standard approach using three-point 
functions. 
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